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Abstract In rapidly rotating stars Coriolis forces and centrifugal deformations mod- 
ify the properties of oscillations; the Coriolis force is important for low-frequency 
modes, while the centrifugal deformation affects mainly p-modes. Here, we dis- 
cuss properties of g- and r-mode oscillations in rotating stars. Predicted frequency 
spectra of high-order g-modes (and r-modes) excited in rapidly rotating stars show 
frequency groupings associated with azimuthal order m. We compare such proper- 
ties with observations in rapidly rotating Be stars, and discuss what is learnt from 
such comparisons. 



1 Oscillations in main-sequence B-stars 

Thanks to OPAL and OP opacity tables fZTl !5|, we now understand that radial 
and nonradial oscillations found in main-sequence B-stars; i.e., j3 Cephei and SPB 
(Slowly pulsating B) stars, are excited by the kappa-mechanism associated with the 
Fe opacity bump at T ^ 2 x lO^K ll28l[37ll23l[20l . Low-order p- and g-modes are ex- 
cited in j3 Cephei stars, while high-order g-modes are excited in SPB stars. Figurefl] 
shows their positions in the HR diagram and predicted instability regions; solid and 
dotted lines are for models without and with core overshooting, respectively. The 
instability regions bounded by solid lines are roughly consistent with j3 Cephei (in- 
verted triangles) and SPB (triangles) stars, which are mostly slow rotators. 

Stellar oscillations give us useful information on the stellar interior which is hard 
to obtain by other means. For j3 Cephei stars having low-order p- and g-modes, 
mode identifications are less ambiguous so that detailed asteroseismic studies are 
possible. Comparing observed frequencies with theoretical ones yields best esti- 
mates of physical parameters as well as the extent of core-overshooting for each 
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Fig. 1 Positions of rapidly rotating Be stars (filled circles) in the HR diagram are shown along 
with the well known B-type main-sequence variables SPB (triangles) and j8 Cep stars (inverted 
triangles). Big filled circles indicate the Be stars with nonradial pulsations detected by space pho- 
tometry from the MOST and CoRoT satellites. Parameters of Be stars are taken mainly from 
1221 1541 . and of SPB stars from 1151 1401 . Also shown are theoretical instability boundaries for p- 
and g-modes. Solid lines for evolutionary tracks and instability boundaries are from models with- 
out core overshooting, while dotted lines come from models with a core overshooting of O.lHp 
(i?p=pressure scale height). Theoretical models were obtained using a standard chemical compo- 
sition of {X,Z) = (0.70,0.02) with OPAL opacity tables 127J. 



Star. In addition, rotational m-splittings of p- and g-modes, which have different 
depth sensitivity, can be used to measure the strength of differential rotation in the 
stellar interior Such asteroseismic analyses have been done for some j3 Cephei stars; 
e.g., V Eri (HEIll^E]), Oph (M^), HD 129929 (US)), j3 CMa (Eg)), 12 Lac 
(HH), and 5 Cet ([2 1). Although results are still somewhat controversial, they seem 
to indicate that the extent of core overshooting in j3 Cephei stars ranges around 0.1- 
Q.3Hp, and the core-to-envelope ratios of rotation rates are approximately 3-5. 

In contrast to the cases of j3 Cephei stars, detailed seismic analysis for SPB stars 
is difficult, because g-mode frequencies are densely distributed, affected strongly by 
stellar evolution, and modified significantly even by moderate stellar rotation; these 
properties make mode identifications difficult (see De Cat [|14J for a review on SPB 
stars). 

Although a detailed frequency fitting between theory and observation for each 
SPB star might be difficult, a collective property as simple as the distribution of SPB 
stars in the HR diagram gives us useful information on the extent of convective-core 
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overshooting. The "red" (or cooler) boundary of the SPB (g-mode) instabihty region 
corresponds to the disappearance of convective cores. The Brunt-Vaisala frequency 
in a radiative dense core is very high and hence the wavelengths of a g-mode become 
very short there, which, in turn, cause strong radiative damping |fT9l l23l l20l . This is 
the reason why the cool boundary from models with a core overshooting of Q.lHp 
(dotted lines in Fig.lTJ is redder than that from models without core overshooting 
(solid lines). Figure[T] shows that the observed SPB stars lay within the instability 
boundary obtained from models without core overshooting. This indicates that, in 
contrast to the cases of j3 Cephei stars, mixing by core overshooting should be weak 
in SPB stars; i.e., in the slowly rotating main-sequence stars of S-SMq. We need 
further accumulation of fundamental parameter data of SPB stars to confirm the 
property. 

Also plotted in Fig.[T] are the positions of rapidly rotating Be stars (filled cir- 
cles). Be stars are B-type stars which have (or had sometime before) Balmer lines 
in emission. The emission lines arise from a circumstellar disk ejected from the star 
due to rapid rotation. Many Be stars are known to show spectroscopic and photo- 
metric variations on various time-scales; in particular, short-term (order of a day) 
photometric (e.g., Balona |7|) and line-profile (e.g., Rivinius & Baade t441 ) vari- 
ations are thought to be caused by radial and nonradial pulsations of Be stars, al- 
though rotational modulations were also suggested for the photometric variations. 
The Be stars with short-term photometric variations are sometimes called A Eri vari- 
ables (Q). Interestingly, Be stars tend to be located in the j3 Cephei (p-mode) and 
SPB (g-mode) instability regions, supporting pulsation (oscillation) origin of the 
short-term variations. Recently, the MOST and CoRoT satellites have found multi- 
periodic light variations in several Be stars, of which amplitude diagrams are shown 
in Fig. 12] The muliti-periodicity strongly supports the explanation that the short-term 
variations of Be stars are caused by pulsations rather than rotational modulations. 

The frequency spectra in the top panels of Fig.|2]for the hot Be stars, Oph and 
HD 49330, show wide ranges of periodicities ranging from about a day to a few 
hours, while the lower panels show relatively long and grouped periods. The dif- 
ference can be understood by the fact that the hottest two stars lay in the j3 Cephei 
instability range (Fig.!!) and p-modes and low order g-modes are excited, while the 
other relatively cooler Be stars lay in the SPB instability range where long-period 
g-modes are excited. Figurefllshows that some cool Be stars including some whose 
variability are detected by CoRoT (large filled circles) lay outside of the SPB in- 
stability range obtained from models without convective core overshooting. This 
indicates that an extensive mixing around the convective core should be occurring, 
if nonradial g-modes are excited in those Be stars. The analyses for HD 181231 
and HD 175869 by Neiner et al.|39| do indicate the presence of such mixing. Fu- 
ture observation of other cool Be stars should tell whether such extensive mixing is 
ubiquitous in rapidly rotating Be stars. 

In contrast to the broad frequency spectra of the hottest Be stars (top panels of 
Fig.|2]i without other remarkable features, the spectra of the relatively cooler Be stars 
show conspicuous frequency groups which are regularly separated. This property 
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Fig. 2 Frequency-amplitude diagrams for liglit variations detected in Be stars by the MOST (left 
panel) and the CoRoT (right panel) satellites. In addition to these stars, MOST detected in the 
late-type Be star j8 CMi small amplitude oscillations with frequencies of 3.26 and 3.28 c/d |46|. 
Frequencies and amplitudes adopted from |52 |-f Oph, 1 10|-HD 127756 & HD 217543, | 531-HD 
163868, |26|-HD49330, |38|-HD 181231, 1 17 |-HD 50209, and |25|-HD 175869. Note that the 
stars obsei-ved by CoRoT tend to be more evolved (with the spectral luminosity class IV) than the 
stars observed by MOST. 



comes from strong rotation effects on high-order g-modes whose intrinsic frequen- 
cies are comparable with or less than the rotation frequencies. 



2 Oscillations in rotating stars 

Stellar rotation affects the oscillations in two ways; through the centrifugal force 
which deforms the equilibrium structure from spherical symmetry, and the Coriolis 
force which represents the angular momentum conservation when matter moves. 
In this section we discuss rotation effects mainly on low-frequency modes (for p- 
modes see e.g., Goupil ll24l '). 



2.1 Coriolis and centrifugal force effects 

The effects of Coriolis and centrifugal forces can be seen from the equation of mo- 
tion, which may be written in the inertial frame as 

d\ 1 

_+vVv = Vp-Vi/A, (1) 

dt p 
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where v is fluid velocity, p matter density, p pressure, and 1/ gravitational poten- 
tial. The velocity consists of oscillation velocity v' and rotation velocity r sin 0X2 e,;,, 
where e,;, is the unit vector of ^ direction. In equilibrium state without oscillations. 
Equation ^ is reduced to 

-rsin0i22en, = -— Vpo-Vv/o, (2) 
Po 

where 

= sin0er + cos0ee = e,|, X Cj, (3) 

is the outwardly pointing unit vector perpendicular to the rotation axis. 

Writing variables in Equation ([T) as a sum of the equilibrium value and (Eule- 
rian) perturbation due to oscillations; e.g., p = po+ p', subtracting Equation (|2]l and 
disregarding non-linear terms with respect to the perturbed quantities, we obtain 

^ + n4l^' + v' ■V{rsmenes) ^ --Vp' + ^Vp ~V\i/ , (4) 
at a(p p 

where the subscript for equilibrium quantities is dropped for ease of notation. 

The relation between oscillation velocity and the Lagrangian displacement ^ is 
somewhat different in a rotating star. Since d£, /At corresponds to the Lagrangian 
velocity of oscillation, the Eulerian oscillation velocity v' is written as 

v' = ^ - ^ V(rsin0i2e^) = ^+^^-^- V(rsin 0i2e^). (5) 

We express the temporal and azimuthal dependence of an oscillation mode as 

ei(«+m0)^ (6) 

where O is the oscillation frequency observed in the inertial frame and m is the az- 
imuthal order of the oscillation. Equation (|6) means that we adopt the convention 
that prograde modes correspond to m < 0; i.e., an oscillation propagates in the pos- 
itive direction of (p (or in the direction of rotation velocity) if m < 0. Then, Eq. (|5]l 
becomes 

v' = iffl<^ -iarsinee^^ • Vi2, (7) 

where 

(0 = a + mn (8) 

is the frequency in the frame rotating with an angular frequency of £2. Note that in 
the case of uniform rotation, co is constant in the stellar interior and can be adopted 
as eigenfrequency. 

Substituting this expression of v' into Equation ([U, we obtain 

1 o' 

- (0^^ +2icoQ{e. X ^) + rsin0(^ • yQ^)es, = V;/ + ^-^Vp - Vi/, (9) 

p p^ 
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where e^. is the unit vector parallel to the rotation axis. The second term in the left 
hand side, 2ia)i2(e,- x corresponds to the Coriolis force. The centrifugal-force 
effect is hidden in Vp, which is, in the non-rotating case, equal to —trGM,p/r^. 
In the presence of rotation the centrifugal force modifies p^'Vp by the order of 
rsinOQ^em and hence breaks spherical symmetry. Thus, the importance of the 
centrifugal force is measured as /cen = r^sm9£2^/GMr, which indicates that the 
effect of centrifugal forces is largest at the equator on the stellar surface where 
/cen = Rlqi2^/GM < 1 with Req being the equatorial radius, and minimum near the 
center where /cen — > ^ sin 9Q^/Gpc <C 1 with Pc being the central density. 
On the other hand, the importance of the Coriolis force can be measured by 

in 

/cor = — , (10) 
CO 

which determines the relative importance between the first two terms on the left 
hand side of equation (|9]i; when /cor > 1 the Coriolis term is larger than the accel- 
eration term (first term) and vice versa. Obviously, for a given rotation frequency 
g-modes are affected more strongly by Coriolis force than p-modes because fre- 
quencies of p-modes in the co-rotating frame are larger than those of g-modes. 

The p-mode frequencies are bounded as (0^ > IQGM/R^ (see e.g.. Cox fTT)), 
while the rotation frequency is limited by < GM/R^. Therefore, for p-modes 
/cor < 1; i-C-, the effects of Coriolis forces are always small for p-modes, while 
the effect of centrifugal force (deformation) can be significant for p-modes because 
the amplitudes of p-modes are confined to outer layers. Actually, Reese et al. B3l 
have shown that p-mode properties are significantly modified in rapidly rotating 
deformed stars. 

For g-modes, on the other hand, the opposite is true; /cor can be larger than unity 
even in a slowly rotating star for high-order g-modes, while the effect of centrifugal 
force is small even in rapidly rotating stars because the amplitude of g-modes is 
confined to inner layers where the centrifugal deformation is small. The insensitivity 
of g-mode frequencies to the centrifugal deformation has been shown numerically 
by Ballot et al. |61 . 

As seen in Equation (|9]), in the absence of rotation the horizontal displacement 
is proportional to ^h{p' / P + V) with 

d Id 
^h = ee— +e^^— — , (11) 
ad ^ sm0 df 

which makes the angular dependence of an oscillation mode representable by a sin- 
gle spherical harmonic Y"''{9,^) (see. Cox 111 11 . Unno et al. |51 1, or Aerts et al. H] 
for details). This simple property is lost in the presence of rotation even if the cen- 
trifugal deformations are neglected, because of the presence of the Coriolis force. 
Then, we express the angular dependence of oscillation by using a sum of terms 
associated with spherical harmonics as 
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^ = '^\sX;^, + HjV^Y;j + T^h^Y;Axer] and p'^tp'^^l", (12) 

where I j = \m\ + 2{j — I) + 1 and Ij — Ij + I ~2I with 7 = for even modes and 
7=1 for odd modes, and J means a truncation length. Other scalar variables are 
expressed in a way similar to p'. The eigenfunction for a scalar variable of an even 
(odd) mode which consists of terms proportional to F," with even (odd) values of 
(Ij — \m\) is symmetric (anti-symmetric) about the equatorial plane. 

As seen in the above equations, we can still assign azimuthal degree m and the 
(even or odd) parity of a mode, but no longer a latitudinal degree. To identify the 
latitudinal behavior of a g-mode we sometimes use in this paper the notation £q; the 
latitudinal degree the mode would have when Q ^Q. When we discuss the obser- 
vational properties, we use the effective degree £ which represents the value of Ij for 
the component with the maximum amplitude on the stellar surface among the other 
components in Eq. (fTSl i. and we use £' for the corresponding toroidal component. 

Note that toroidal components are needed in representing the displacement vec- 
tors, because the toroidal velocity fields couple, through the Coriolis force, with 
spheroidal velocity fields which generate variations in density, pressure, tempera- 
ture, etc. 



2.2 r-modes 

Purely toroidal motions in a non-rotating spherical star do not disturb stellar struc- 
ture so that no restoring force nor oscillations arise. In the presence of rotation, 
however, the latitudinal gradient of the vertical component of the angular frequency 
of rotation i2 cos 9 causes a restoring force for toroidal motions so that toroidal os- 
cillations occur; we call these oscillations r-modes (or global Rossby modes). All the 
r-modes are retrograde in the co-rotating frame (see e.g., Pedlosky [41], Saio |45|). 
If we assume that the displacement is purely toroidal and neglect the centrifugal 
deformation, the toroidal oscillation would have the limiting frequency, Oriim given 
as 

'^■'™ = f(FTT)' ^^^^ 

in the co-rotating frame. Actual frequencies of r-modes in the co-rotating frame 
deviate from Oriim; the deviation is larger for higher radial order modes (Provost 
et al. ||42l. Saio B31 ). The deviation arises because the Coriolis forces associated 
with toroidal motion generate spheroidal motions; i.e., vertical motion and density 
perturbations are generated, and hence buoyancy plays a role in the restoring forces 
for r-modes. Because of these effects, r-modes are sometimes called "mixed modes" 
(Townsend Il49l ) or "q-modes" (quasi g-modes; Savonije 147 1). The coupling with 
spheroidal motions makes it possible for r-modes (mainly toroidal oscillations) to be 
excited thermally by the kappa-mechanism in the same way as g-modes are excited. 
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Stability analyses by Townsend ||49l . Savonije ||47l . and Lee ||30l . indicate that some 
r-modes are actually excited by the kappa-mechanism in intermediate mass (around 
3-8M0) main sequence stars, in which g-modes are also excited (SPB stars). 



2.3 Latitudinal amplitude distributions 

Latitudinal distributions of oscillation amplitudes are expressed as shown in Equa- 
tion ([T2]i, which indicates a dependency on the ratio /cor- The distributions are im- 
portant in understanding the visibilities of excited modes. 

Figure[3] shows latitudinal distributions of the amplitude of p' as a function of 
COS0 for selected g-modes of m = ±1 and £0 = I, 2, in a 4.5Mq ZAMS model 
of rigid rotation at an angular frequency of .220 GM / (or with a period of 
0.96 days, Veq — 130km/s). Although the adopted rotation rate is moderate, the Cori- 
olis force significantly affects the amplitude distribution of g-modes, because fre- 
quencies are so low that the parameter /cor = 2Q /(O is larger than around 1. These 
eigenfunctions have been obtained by the method of Lee & Baraffe \ 33\ without us- 
ing the 'traditional' approximation in which the horizontal component of rotational 
angular velocity, Q sin 9eg is neglected. The results are very similar to those ob- 
tained by using the traditional approximation (see, e.g., Lee & Saio f34l), indicating 
that the approximation is well suited for low-frequency oscillations in which the 
horizontal velocity is much larger than the radial one. 




o.a 0.4 0.6 0.8 1 0.2 0.4 0.6 o.a 



cos 6 cos 6 

Fig. 3 Amplitude distributions for scalar variables of g-modes with m = ±1 and f.o = I and 2 as 
a function of cos 9, where 9 is co-latitude. Left and right panels are for prograde (m = —1) and 
retrograde (m = 1) g-modes, respectively. Solid and dashed lines are for even {('q = 1) and odd 
(^0 = 2) modes, respectively. A number along each line indicates the value of /cor = 2X2 /a> for the 
mode. Note that an additional nodal line appears for retrograde modes when /cor > 1, which would 
reduce the visibility. Dotted line shows the associated Legendre function (cos 6). 
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O.a 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 



COS 6 COS 

Fig. 4 The same as Fig.[3]but for g-modes of m = (left panel) and r-modes of m = 1 (right panel). 
In the left panel, solid and dashed lines are = 2 (even) and = 1 (odd) modes, respectively. In 
the right panel, solid lines are for even modes with (. = (' — 1 = 1, and dashed lines are for odd 
modes wixh t = (' + \ =2. Dotted lines shows the Legendre function of Pj (cos 9 ) (left panel) and 
(cos 6 ) (right panel). 



Generally, amplitudes of oscillation modes with large /cor tend to be confined to 
equatorial regions. This property is consistent with the behavior of the eigenvalue 
X of the Laplace tidal equation (e.g., Lee & Saio |34l and Aerts et al. |T1), which 
governs angular dependence of oscillation in the traditional approximation. As /cor 
increases, A increases rapidly except for sectoral prograde modes (^o = — m); a large 
X corresponds to a large effective degree £, i.e., X £o(^o + 1) as Q — > 0. 

Solid lines in Fig. [3] are for = \m\ = 1; i.e., sectoral modes, which have no 
latitudinal nodal lines when X2 = 0. The prograde modes (m = — 1) keep this prop- 
erty even for a large /;or, while for retrograde modes {m = 1) a latitudinal nodal 
line appears for the cases with /cor > 1, which reduces visibility of the mode due to 
cancellation on the surface. 

The left panel of Fig.|4]shows latitudinal amplitude distributions of axisymmetric 
(m — 0) g-modes (dashed lines for = 1 modes, and solid lines for = 2 modes). 
Similarly to non-axisymmetric modes, amplitudes of axisymmetric modes tend to 
be confined in equatorial zones as /cor increases; in particular for = 1 modes the 
amplitude weight shifts from the pole to low-latitudes. 

The right panel of Fig. |4] shows amplitude {p') distributions of selected r-modes 
of m = 1 in the same model as in Fig. [3] Red solid lines are for even r-modes with 
£' = £-1-1=2 and blue dashed lines are for odd r-modes with £' = £—1 = 1, where 
£' and £ are the latitudinal degrees for the dominant toroidal component and for the 
corresponding spheroidal component, respectively. (The parity refers to the property 
of scalar variables.) 

For r-modes of m = 1 (right panel of Fig.|4]l, a low-order mode close to the limit- 
ing frequency Oriim (Eq. ( fTsT i) and a relatively high-order mode are shown for each 
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parity. (For low-order r-modes of m — 1, f^oY ^ + I)-) The amplitude distri- 
bution of a higher order r-mode tend to be more confined to the equatorial region, 
which is similar to high-order g-modes having large /cor- An important difference 
from retrograde g-modes is that no additional latitudinal nodal-line appears for r- 
modes. The stability analyses for r-modes by Townsend 1491 and Lee ||30) indi- 
cate that odd r-modes are more easily excited in B-type main-sequence stars, while 
Savonije [47 1 found even (symmetric with respect to the equator) r-modes to be ex- 
cited in some models. We expect odd r-modes to be detected unless the inclination 
angle between rotation axis and the line-of-sight is close to 90°. 



2.4 Expected frequency ranges of g- and r-modes 

High order g-modes (Gautschy & Saio [1231 . Dziembowski et al. fSOl) and r-modes 
(Townsend |49|, Savonije pTl . Lee ll30l ) are excited by the kappa-mechanism at 
the Fe-opacity bump in intermediate mass (around main-sequence stars. 

Some g-modes are, however, damped in rapidly rotating stars (/cor > 1) due to mode 
couplings. The mode coupling occurs when two normal modes with the same m 
and parity but with different io have similar oscillation frequencies, as these two 
modes are no longer independent in the presence of the Coriolis force. Lee |29| 
has found that the mode couplings tend to damp retrograde g-modes. The effects 
of the mode couplings are discussed in detail in Aprilia et al. [3 1 and Lee lf3n[32l . 
The damping effects can also be significant on prograde g-modes except for sectoral 
modes {£o = —m). Although some prograde tesseral (£o > —m) g-modes might still 
remain excited in rapidly rotating stars, the effective degrees (£) of these modes 
might be too high to be detectable. Therefore, most visible modes should be low- 
degree prograde sectoral (£0 = —m) g-modes and (retrograde) r-modes in rapidly 
rotating intermediate mass stars. 

For high order g-modes in rapidly oscillating stars o ^ \m\Q with o being 
frequencies in the co-rotating frame so that we expect to observe these modes at 
a = \(0 — mQ.\ ~ |m|X2. 

Since sectoral prograde g-modes are expected most visible as discussed above, 
we expect groups of frequencies slightly above Q. and 2Q. for m= —1 and —2 
modes, respectively. (Modes with higher \m\ are expected to be less visible due to 
cancellation.) 

For r-modes with m — I only odd modes (£' = 1) seem to be excited in B-type 
stars, they have co ^ Q, whose observational frequencies are very small <T ^ X2. 
Symmetric (even) r-modes of ;« = 2 (^ = 2) are also excited in some cases, for 
which £' = 3, hence cOcUm = ^^2. Therefore, frequencies in the observer's frame are 
slightly larger than | Q . 

In summary, we expect to observe groups of frequencies at approximately 



g-modes: X2, 2X2, and r-modes: 0, -Q, 



(14) 
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Fig. 5 Growth rates and azimuthal order m of low-degree (£ < 2) modes excited in 4.5MiV, main- 
sequence models with various rotation frequencies, where the results were obtained by the method 
of Lee & Baraffe ||33i|. The cyclic rotation frequency is indicated by a vertical dash-dotted line at 
each panel, and Q is the corresponding angular frequency of rotation normalized by ^JCM/W 
with R being the mean stellar radius; equatorial rotation velocity is also shown. Horizontal axis 
indicates frequency in the observer's frame. 



where we have assumed arbitrarily that modes with effective degrees I <2 are visi- 
ble. 

Figure|5] shows growth rates and azimuthal order m of the low degree (£ < 2) 
modes excited in a A.5Mq main-sequence model of rigid rotation at various speeds. 
The abscissa indicates frequencies in the observer's frame. Red (solid) and blue 
(dotted) lines are for even and odd modes, respectively. In the model of the lowest 
rotation frequency (12 = 0.073, bottom panel), no r-modes are excited, and we see 
no appreciable effects of rotation on the stability of g-modes; rotation only disperses 
frequencies by the effect of transformation from co-rotating to inertial frames; a = 
0) — mi2. 

In the second lowest rotation case with a normalized angular frequency of rota- 
tion Q. — 0.15, m = 1 odd r-modes with ^' = ^ — 1 = 1 are excited; they have very 
small frequencies in the observer's frame, because (Oram = ^ (Eq. (fTSl )). 

As the rotation frequency increases (Q > 0.2), retrograde (m > 0) g-modes tend 
to be damped, while prograde g-modes remain excited. Due to the damping of ret- 



12 



Hideyuki Saio 



rograde g-modes and the increasing effect of —mQ, frequency groupings become 
conspicuous as the rotation frequency increases. In a sufficiently rapidly rotating 
case (as in the top panel), well populated frequency groups of prograde sectoral 
(£ = —m) g-modes are formed around 1.2Q and 2.3Q corresponding to m = — 1 
and —2, respectively. 

The r- modes also form groups, because the frequency deviation \cOr\im — ft) | of 
an r-mode is usually much smaller than mQ. The group at smallest frequencies 
corresponds to m = 1 anti-symmetric r-modes £' = £—! ~ 1 . In the fastest rotation 
case (top panel) all the excited retrograde (m > 0) modes are r-modes. Note that in 
this model m ~ 2 r-modes with £ — £' — ! =2 symmetric with respect to the equator 
are excited; their observational frequencies are slightly larger than | Oiiim — inQ \ = 
|X2. The excitation of symmetric m = 2 r-modes in a sufficiently rapid rotator agree 
with the results of Savonije ['47). 

Figure|5]shows that frequency groupings appear even in stars whose rotation rates 
are considerably less than the critical rate. This is consistent with the finding of 
Balona et al. |8) from Kepler data that several B stars with relatively large V sin; 
show frequency groupings. 



3 Comparisons with frequency groups of Be stars 

The characteristic frequency groupings obtained for several Be stars (Fig.|2]) were 
fitted with models (e.g., Cameron et al.|10|, Neiner et al.[39|). Such fittings give 
estimates for rotation frequencies and in some cases, the constraint to the internal 
mixing in rapidly rotating stars. In this section, we discuss what we learn from model 
fittings to low-frequency oscillations of Be stars. 



3.1 HD 50209 

We will discuss, as an example, a model fitting to low-frequency oscillations of 
the late-type (B8IVe) Be star HD 50209 detected by CoRoT (| 17|); the ampHtude 
spectrum consists of five frequency groups around O.lc/d, 0.6-0.8c/d, 1.5c/d, 2.2c/d 
and 3c/d (or six groups if the 0.6-0.8c/d group is separated into two groups) as 
shown in Fig.|2]and in the left-bottom panel of Fig.|6] 

The fundamental parameters of HD 50209 have been derived by Diago et al. ifTTl 
as logTeff = 4.134 ±0.051, logL = 3.02 ±0.39, and logg = 3.56 ±0.11 (with 
Q /i2crit = 0.90). The estimated position of HD 50209 in the HRD is shown in the 
right panel of Fig.|6]by a filled square with en^orbars. At the center of the parameter 
ranges the critical rotation rate is around .7c/d. Since prograde g-modes of m — — l 
form a group around frequencies slightly larger than the rotation frequency, the 0.6- 
O.Sc/d group should be fitted by assuming a rotation frequency of approximately 
0.6c/d. The rotation frequency should not be far from the critical rate; i.e., the ra- 
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M= 5.5 log L= 3,240 log Te=_4,132 log R- 0,880 
•^rot =0,605 e/d; 0= 0,631 
B 4 k |m| S 3; 1 'S S ' ' ' ' ' ' ' ' ' ^ 




Frequency(c/d) ^.^ ^ 

Fig. 6 Left panel: Frequencies of excited modes in a model of a S.SMt, star rotating at a rate 
of 0.605 c/d are compared with observed frequencies of HD 50209 (bottom panel) adopted from 
Diago et al. 1 17 1; the rotation frequency of the model is indicated by the dash-dotted vertical line. 
The top panel shows azimuthal order m (prograde modes con'espond to m < 0) and the middle 
panel shows growth rates of excited modes. Solid (red) and dashed (blue) lines are even and odd 
modes respectively. Right panel: The positions of HD 50209 (filled square with error bars) on the 
HRD and two models (filled circles) that are reasonably consistent with the observed frequency 
groups; a fit by one of the models is shown in the left panel. Evolutionary tracks are calculated 
with a standard composition of {X,Z) = (0.70,0.02). The thick two tracks for 5.0Mq and 5.5Mq 
are from evolution including overshooting of 0.35//,, around the convective core; the other tracks 
are without overshooting. 

dius should be sufficiently large, otherwise no clear frequency groupings would be 
expected (see Fig.|5]l. Two models of S.OM© and 5.5M0 stars rotating at a rate of 
0.6 c/d (0.007 mHz) were examined; their positions on the HRD are shown by filled 
circles in Fig.|6] In both cases a substantial core overshooting (more than Q.3Hp) 
must be included to meet the above requirements for T^ff and luminosity within the 
spectroscopic estimates. The requirement of substantial core overshooting for the 
model of HD 50209 is common to the models for other CoRoT Be stars, HD 175869 
and HD 181231 as discussed in Neiner, et al. |[39l . 

Both models produce similarly good fits to observed frequency groups of HD 50209. 
A comparison with the 5.5Mp model is shown in the left panel of Fig.|6] where the 
method of Lee & Baraffe |33 | was used for the stability analysis. To fit the three 
(0.6-0.8, around 1.5 and 2.2 c/d) groups with g-modes, we have to consider az- 
imuthal orders up to 3 (\m\ < 3). (The frequencies at around 3 c/d are probably 
harmonics of large amplitude frequencies at approximately 1 .5 c/d as indicated by 
Diago et al yjj. ) 
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One can interpret the smallest frequency group at around . 1 c/d as due to m = 
1 anti-symmetric r-modes. Diago et al. ifTTl estimate an inclination angle of / ^ 
60° for HD 50209 for a rotation frequency of 90-95% of the critical rate. With 
this inclination angle, r-modes would suffer from no serious surface cancellation 
(see right panel of Fig.|4|i, and would be detectable if the r-modes produce enough 
temperature variations. It is interesting to note that a very low-frequency group is 
common in SPBe stars that show g-mode pulsations, but it does not appear in hot Be 
stars Oph and HD 49330) with p-mode pulsations (Fig.|2]i. This agrees with the 
theoretical prediction for the excitation range of r-modes on the HRD as obtained 
by Townsend |49|, and supports the identification of the very low frequency groups 
in SPBe stars as r-modes. In this model, symmetric r-modes of m = 2 and 3 are 
excited, but they have no observational correspondence. 



3.2 Rotation rates of Be stars 

Figure|7]plots rotation frequencies obtained by fitting theoretical predictions for fre- 
quency groups against those observed for several Be stars observed by MOST and 
CoRoT (a similar diagram is shown in |[TOl ). The lower panel shows the rotation 
frequency of the Be stars (with error bars) as a function of the effective temperature. 
The Be stars observed by CoRoT tend to have low rotation frequencies compared 
to those observed by MOST. This is due to the fact that CoRoT-observed Be stars 
tend to be more evolved, having spectral luminosity class IV. This apparent differ- 
ence disappears if the rotation frequency is normalized by ^ GM/R^ as plotted in 
the top panel of Fig.|2] where R is mean radius. With this normalization critical ro- 
tation corresponds to GM/R^ 0.7. This indicates that all Be stars shown 
can be considered rotating nearly critically, which agrees with recent spectroscopic 
estimates by Townsend et al. |50| and Cranmer 1 12|. 

Dashed and dot-dashed lines in the bottom panel of Fig. |7] show critical rotation 
frequency as a function log Tf.ff for ZAMS and TAMS (terminal age main-sequence) 
models, respectively. For TAMS models, critical rotations with and without substan- 
tial core overshooting cases are shown. Obviously, substantial overshooting would 
need to be assumed for models to be consistent with the CoRoT observations of Be 
stars. 

Also plotted in this figure are rotation frequencies of X Eh variables (filled 
squares); a group of Be stars which show short-term (order of a day) periodic light 
variations known from ground-based observations (e.g., Balona |7|). Folded light 
curves of A Eri variables, some of which show double-wave character, are consis- 
tent with the grouped frequency distributions obtained for relatively late-type Be 
stars monitored by MOST and CoRoT satellites. In fact, the folded light curve of 
HD 181231 from CoRoT observations (Neiner et al. |38 1), for example, is similar to 
typical light curves of X Eh variables. Since a model with a rotation frequency about 
20% smaller than the frequency at the group of m = — 1 modes fits well to such a 
characteristic frequency distribution, rotation frequencies of A Eri variables are as- 
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Fig. 7 Effective temperature versus rotation frequencies estimated by model fittings for Be stars 
observed by MOST and CoRoT (large circles with enm bars). The top panel shows rotation fre- 
quency normalized by GM/E?, where R is the mean radius. With this normalization the critical 
rotation frequency corresponds to approximately 0.7. Filled squares are X Eri variables whose ro- 
tation frequencies are assumed to be 80% of the observed frequencies listed in Balona |7|. The 
dashed line indicates the critical rotation for zero-age main-sequence models as a function of 
log Teff . Dotted-dashed lines show critical rotation of TAMS (terminal-age main sequence) models 
without and with substantial convective core overshooting. 



sumed as 80% of the observed frequencies, and those values are plotted in Fig.|7] 
This figure shows that the distribution of A Eri rotation rates are consistent with 
Be star models rotating nearly critically if a substantial core-overshooting of around 
035Hp is assumed to occur in these stars. Since models without core overshooting 
are consistent with the slowly rotating SPB stars, rapid rotation seems the cause of 
the substantial mixing needed around the convective core (Neiner et al. i39i ). 



4 Conclusion 



Rapid stellar rotation modifies properties of oscillations in complex ways and makes 
comparisons between observational and theoretical results difficult. In particular. 
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low order p- and g-modes excited in early type Be stars in the j3 Cephei instability 
region show complex frequency spectra affected considerably both by deformation 
of the equilibrium structure and by the Coriolis force. The frequency distributions of 
these stars are such that it would seem impossible to identify modes and so compare 
with theoretical results. 

In contrast to the p-modes, rapid rotation helps us to identify the azimuthal de- 
gree m of low-frequency modes, because observational frequencies of these modes 
form groups; i.e., g-modes with the same m have similar frequencies. This prop- 
erty can be used to identify m and to estimate the rotation frequency of each star. 
The rotation rates thus obtained for several Be stars indicate that Be stars rotate 
nearly critically, which confirms conclusions drawn from spectroscopic analyses. 
In addition, in order to fit observations of some Be stars, models with substantial 
core-overshooting are needed. This indicates that extensive internal mixing exterior 
to the convective core occurs in rapidly rotating stars as discussed in Neiner et al. 
||39ll . Furthermore, it might be possible in the near future to compare each g-mode 
(and r-mode) frequency (or period spacings) with theory in order to obtain detailed 
information on the stellar interior Rapid rotations are advantageous in this case be- 
cause we can identify the azimuthal order m while effects of centrifugal deformation 
remain small, as shown by Ballot et al. [6i. 

Acknowledgements I am grateful to Alfred Gautschy and Umin Lee for helpful comments. 



References 

1. Aerts, C, Christensen-Dalsgaard, J., Kurtz, D.W.: Asteroseismology. Springer, Heiderberg 
(2010) 

2. Aerts, C, Marchenko, S.V., Matthews, I.M., Kuschnig, R., et al: S Ceti is not monoperiodic: 
Seismic modeling of a /3 Cephei star from MOST space-based photometry. ApJ, 642, 470- 
477 (2006) 

3. Aprillia, Lee, U., Saio, H.: Stability of g modes in rotating B-type stars. MNRAS, 412, 2265- 
2276 (2011) 

4. Ausseloos, M., Scuflaire, R., Thoul, A., Aerts, C: Asteroseismology of the /3 Cephei star v 
Eridani: massive exploration of standard and non-standard stellar models to fit the oscillation 
data. MNRAS, 355, 352-358 (2004) 

5. Badnell, N.R., Bautista, M.A., Bufler, K, Delahaye, R, Mendoza, C, Palmeri, P., Zeippen, 
C.J., Seaton, M.I.: Updated opacities from the Opacity Project. MNRAS, 360, 458^64 
(2005) 

6. Ballot, 1., Lignieres, Prat, V., Rees, D. R., Rieutord, M. : 2D computations of g modes. In: 
Shibahashi, H. (ed.) Progress in solar/stellar physics with helio- and asteroseismology, in 
press (2012) arXiv: 1 109.6856v2[astro-ph.SR] 

7. Balona, L.A.: Test of the pulsation and starspot models for the periodic Be stars. MNRAS, 
277, 1547-1554 (1995) 

8. Balona, L.A., Pigulski, A., De Cat, P., Handler, G., et al: Kepler observations of the vaiiability 
in B-type stars. MNRAS, 413, 2403-2420 (2011) 

9. Briquet, M., Morel, T. Thoul, A., Scuflaire, R., Miglio, A., Montalban, J., Dupret, M.-A., 
Aerts, C: An asteroseismic study of the /3 Cephei star 6 Ophiuchi: constraints on global 
stellar parameters and core overshooting. MNRAS, 381, 1482-1488 (2007) 



Prospects for asteroseismology of rapidly rotating B-type stars 



17 



10. Cameron,C., Saio, H., Kuschnig, R., Walker, G.A.H.: MOST detects SPBe pulsations in HD 
127756 and HD 217543: Asteroseismic rotation rates independent of vsini. ApJ, 685, 485- 
507 (2008) 

11. Cox, J.P.: Theory of Stellar Pulsation. Princeton Univ.Press, Princeton (1980) 

12. Cranmer, S.R.: A statistical study of threshold rotation rates for the formation of disks around 
Be stars. ApJ. 634, 585-601 (2005) 

13. Daszyriska-Daszkiewicz, J., Walczak, P.: Complex asteroseismology of the /3 Cep/slowly pul- 
sating B-type pulsator v Eridani: constraints on opacities. MNRAS, 403, 496-504 (2010) 

14. De Cat, P.: Observational asteroseismology of slowly pulsating B stars. Comm.Astero., 150, 
167-174 (2007) 

15. De Cat, P., Aerts, C: A study fo bright southern slowly pulsating B stars. II. The intrinsic 
frequencies. A&A, 393, 965-981 (2002) 

16. Desmet, M., Briquet, M., Thoul, A., Zima, W., et ah: An asteroseismic study of the P Cephei 
star 12 Lacertae: multisite spectroscopic observations, mode identification and seismic mod- 
elling. MNRAS, 396, 1460-1472 (2009) 

17. Diago, P.D., Gutierrez-Soto, J., Auvergne, M., Fabregat, J., et al.: Pulsations in the late-type 
Be star HD 50209 detected by CoRoT. A&A, 506, 125-131 (2009) 

18. Dupret, M.-A., Thoul, A., Scuflaire, R., Daszyriska-Daszkiewicz, J., Aerts, C, Bourge, P.- 
O., Waelkens, C, Noels, A.: Asteroseismology of the /3 Cep star HD 129929 11. Seismic 
constraints on core overshooting, internal rotation and stellar parameters. A&A, 415, 251- 
257 (2004) 

19. Dziembowski, W.A.: Nonradial Oscillations of Evolved Stars I. Quasiadiabatic Approxima- 
tion. AcA, 21, 289-306 (1971) 

20. Dziembowski, W.A., Moskalik, P., Pamyatnykh, A. A.: The opacity mechanism in B-Type 
Stars - 11. Excitation of high-order g-modes in main sequence stars. MNRAS, 265, 588-600 
(1993) 

21. Dziembowski, W.A., Pamyatnykh, A.A.: The two hybrid B-type pulsators: v Eridani and 12 
Lacertae. MNRAS, 385, 2061-2068 (2008) 

22. Fremat, Y., Neiner, C, Hubert, A.-M., Floquet, M., Zorec, J.,Janot-Pacheco, E., Renan de 
Medeiros, J.: Fundamental parameters of Be stars located in the seismology fields of COROT. 
A&A, 451, 1053-1063 (2006) 

23. Gautstchy, A., Saio, H.: On non-radial oscillations of B-type stars. MNRAS, 262, 213-219 
(1993) 

24. Goupil, M.J.: Effects of Rotation on Stellar p-Mode Frequencies. Lect. Notes Phys. 765, 45- 
99 (2009) 

25. Gutierrez-Soto, J., Floquet, M., Samadi, R., Neiner, C, et al.: Low-amplitude variations de- 
tected by CoRoT in the B811e star HD 175869. A&A , 506, 133-141 (2009) 

26. Huat, A.-L., Hubert, A.-M., Baudin, R, Floquet, M. et al.: The BO.SlVe CoRoT target HD 
49330 1. Photometric analysis from CoRoT data. A&A, 506, 95-101 (2009) 

27. Iglesias, C.A., Rogers, RJ.: Updated Opal Opacities. ApJ, 464, 943-953 (1996) 

28. Kiriakidis, M., El Eid, M.F., Glatzel, W.: Heavy element opacities and the pulsations of j8 
Cepheid stars. MNRAS, 255, lp-5p (1992) 

29. Lee, U.: Pulsational stability of g-modes in slowly pulsating B stars. ApJ, 557, 311-319 
(2006) 

30. Lee, U.: r modes of slowly pulsating B stars. MNRAS, 365, 677-687 (2006) 

31. Lee, U.: Pulsation in rapidly rotating stars. Comm.Astero, 157, 203-208 (2008) 

32. Lee, U.: Low frequency oscillations in rotating stars, in these proceedings (2012) 

33. Lee, U., Baraffe, 1.: Pulsational stability of rotating main sequence stars: the second order 
effects of rotation on the nonadiabatic oscillations. A&A, 301, 419^32 (1995) 

34. Lee, U., Saio, H.: Low-frequency nonradial oscillations in rotating stars. I.Angular depen- 
dence. ApJ, 491, 839-845 (1997) 

35. Lovekin, C.C., Goupil, M.-J.: Rotation and convective core overshoot in Q Ophiuchi. A&A, 
515, A58 (2010) 

36. Mazumdar, M., Briquet, M., Desmet, M., Aerts, C: An asteroseismic study of the jS Cephei 
star jS Canis Majoris. A&A, 459, 589-596 (2006) 



18 



Hideyuki Saio 



37. Moskalik,P., Dziembowski, W.A.: New opacities and the origin of the Beta Cephei pulsation. 
A&A, 256, L5-L8 (1992) 

38. Neiner, C, Gutierrez-Soto, J., Baudin, R, de Batz, B.: The pulsations of the B5IVe star HD 
181231 observed with CoRoT and ground-based spectroscopy. A&A, 506, 143-151 (2009) 

39. Neiner, C, Mathis, S., Saio, H., Lovekin, C, Eggenberger, P., Lee, U.: Seismic modelling of 
the late Be stars HD 181231 and HD 175869 observed with CoRoT: a laboratory for mixing 
processes. A&A, 539, A90 (2012) 

40. Niemczura, E.: The main parameters of SPB stars on the basis of lUE spectra. In: Aerts, C, 
Bedding, T.R., Christensen-Dalsgaard, J. (ed.) Radial and Nonradial Pulsations as Probes of 
Stellar Physics, ASP Conf. Ser. 259, 230-231 (2002) 

41. Pedlosky, J.: Geophysical Huid Dynamics. Springer- Verlag, Berlin (1979) 

42. Provost, J., Berthomieu, G., Rocca, A.: Low frequency oscillations of a slowly rotating star - 
Quasi toroidal modes. A&A, 94, 126-133 (1981) 

43. Reese, D.R., MacGregor, K.B., Jackson, S., Skumanich, A., Metcalfe, T.S.: Pulsation modes 
in rapidly rotating stellar models based on the self-consistent field method. A&A, 506, 189- 
201 (2009) 

44. Rivinius, T, Baade, D., Stefl, S.: Non-radially pulsating Be stars. A&A, 411, 229-247 (2003) 

45. Saio, H.: R-mode oscillations in uniformly rotating stars. ApJ, 256, 717-735 (1982) 

46. Saio, H., Cameron, C, Kuschnig, R., Walker, G.A.H., et al.: MOST detects g-modes in the 
late-type Be star j3 Canis Minoris (B8 Ve). ApJ, 654, 544-550 (2007) 

47. Savonije, G.J.: Unstable quasi g-modes in rotating main-sequence stars. A&A, 443, 557-570 
(2005) 

48. Suarez, J.C., Moya, A., Amado, P.J., Martm-Ruiz, Rodriguez-Lopez, Garrido, R.: Seismology 
of P Cephei stars: Differentially rotating models for interpreting the oscillation spectrum of 
V Eridani. ApJ, 690, 1401-1411 (2009) 

49. Townsend, R.H.D.: Kappa-mechanism excitation of retrograde mixed modes in rotating B- 
type stars. MNRAS, 364, 573-582 (2005) 

50. Townsend, R.H.D., Owocki, S.P., Howarth, LD.: Be-star rotation: how close to critical? MN- 
RAS, 350, 189-195 (2004) 

51. Unno, W., Osaki, Y., Ando, H., Saio, H., Shibahashi, H.: Nonradial pulsations of stars. Univ. 
Tokyo Press, Tokyo (1989) 

52. Walker, G.A.H., Kuschnig, R., Matthews, J.M., Reegen, P. et al.: Pulsations of the Oe star ^ 
Ophiuchi from MOST satellite photometry and ground-based spectroscopy. ApJ, 623, L145- 
L148 (2005) 

53. Walker, G.A.H., Kuschnig, R., Matthews, J.M., Cameron, C. et al.: MOST detects g-modes 
in the Be star HD 163868. ApJ, 635, L77-L80 (2005) 

54. Zorec, J., Fr6mat, Cidale, L.: On the evolutionary status of Be stars I. Field Be stars near the 
Sun. A&A, 441, 235-248 (2005) 



